International Journal for Numerical Methods in Biomedical Engineering

Classical neuronal cable theory relies on quasi-static electric field approximations and neglects magnetic induction, Lorentz force coupling, and transient electromagnetic currents, limiting its ability to fully characterize action potential propagation within geometrically branched axons and dendrites. This work develops a coupled Maxwell-electromagnetic cable framework by integrating finite-difference time-domain (FDTD) solutions of Maxwells equations with extended Hodgkin-Huxley and Fitzhugh-Nagumo membrane dynamics, incorporating magnetic gating perturbations, electromagnetic trans-membrane currents IEM, and nanoscale quantum corrections for thin neural segments. Controlled propagation experiments are designed to quantify deviations from standard cable predictions across asymmetric and symmetric axonal bifurcation geometries. Numerical results demonstrate that inductive magnetic effects lower the critical branch radius for junction conduction failure and break symmetric action potential invasion in geometrically identical child branches under external transverse magnetic fields. An electromagnetic corrected geometric ratio GREM is proposed to revise impedance-matching conditions at branch points, accounting for size-dependent axial current imbalance induced by magnetic and displacement currents. Parent axon conduction velocity deviates substantially from the canonical [Formula] scaling law when electromagnetic feedback and quantum charge distributions are included, triggering early signal blockage at large cable diameters. Collectively, this study establishes that quasi-static cable models underestimate electromagnetic corrections to propagation speed, waveform shape, and bifurcation transmission fidelity; the coupled Maxwell-cable framework provides a comprehensive multi-physics tool for modeling electrodynamic signal behavior in complex neuronal architectures.

Background: Cerebral autoregulation stabilizes cerebral blood flow over a range of cerebral perfusion pressures, but the precise shape of the pressure-flow relationship remains debated. The classical triphasic pressure-flow relationship was recently challenged by experiments demonstrating a quadriphasic response, hypothesized to arise from vessel-size-dependent pressure-diameter responses. We tested this hypothesis and investigated whether these size-dependent responses originate from heterogeneity in smooth muscle cell (SMC) abundance, SMC behavior, or neither. Methods: We developed a computational multi-scale physics-based model of cerebral autoregulation linking SMC activity to vessel-scale diameter regulation and organ-scale blood flow. Four scenarios were evaluated: passive vessels, homogeneous SMC abundance and behavior, heterogeneous SMC abundance, and heterogeneous SMC behavior. Predicted pressure-diameter responses and pressure-flow relationships were compared across scenarios and against experimental observations. Results: In contrast to passive vessels, homogeneous SMC activation produced partial flow stabilization, highlighting the key role of SMCs in autoregulation. However, only heterogeneous SMC behavior reproduced the experimentally observed vessel-size-dependent trends in pressure-diameter responses. This scenario also showed the best agreement with the experimental organ-scale pressure-flow relationship (R-squared = 0.93, nRMSE = 5.96%). Conclusion: The model suggests that vessel-size-dependent SMC behavior underlies vessel-size-dependent pressure-diameter responses and shapes the relationship between cerebral perfusion pressure and cerebral blood flow.

One goal of Scientific Machine Learning (SciML) is to advance traditional scientific computing frameworks with modern machine learning tools. This includes extending established methods, such as the finite element method, with cardiac function applications due to their complexity and need for very rapid execution times for real time clinical use. In this work, we present an advanced form of the Neural Network Finite Element (NNFE) method specialized for cardiac simulations, termed CARDIAX-NNFE. The NNFE method learns the parameter-to-displacement field map by training over the residual of the hyperelastic material PDE, using the domain represented by finite elements. The implementation is developed in Python using JAX to leverage its automatic differentiation, highly parallel GPU, and JIT-compilation capabilities. To demonstrate CARDIAX-NNFE effectiveness, we trained full cardiac pressure-volume responses using a simplified heart model, spanning the entire cardiac physiological functional range. Results indicated the ability to simulate a family of pressure-volume solutions with average nodal positional error of 0.023 mm and maximal error of 0.054 mm, with a single complete PV loop evaluated in 0.002 seconds. The CARDIAX-NNFE software platform thus provides for a robust platform for cardiac functional simulations. Moreover, it provides the structure for residual-based SciML methods, which can apply to a variety of physics-based biomedical problems that require high execution speed for clinical applications.

The temporal resolution of medical imaging sequences used to drive patient-specific computational fluid dynamics (CFD) simulations remains limited, typically providing 10-20 frames per cardiac cycle. Therefore, temporal interpolation to reconstruct left atrial (LA) wall motion and boundary conditions is required, but its impact on hemodynamic predictions has not been systematically characterized. To investigate this, we constructed high-temporal-resolution reference wall-motion data using electromechanical (EM) simulations on five patient-specific atrial geometries with a history of atrial fibrillation. We then generated temporally downsampled datasets to emulate clinical frame rates (5, 10, 20, and 40 frames per cycle) and performed CFD simulations to isolate the effects of temporal undersampling on hemodynamic metrics. The focus was placed on kinetic energy, KE, and residence time, TR, particularly in the left atrial appendage (LAA), where thrombosis is most likely to occur. We employed an immersed boundary method to prescribe the wall motion and computed blood TR through a passive scalar transport equation. Results indicate that while global LA hemodynamic indices were marginally affected by the frame rate (errors < 9%), LAA metrics were more sensitive with errors up to 31% compared to reference values. The results based on 20 and 40 frames per cycle yielded favorable agreement with reference results, while 5-and 10-frame reconstructions showed larger, though not systematically biased, deviations from the reference. Importantly, patient ranking by blood-stasis indices was largely preserved. The analysis suggests that patientspecific LA reconstructions derived from dynamic CT imaging provide a reliable basis for estimating LAA blood-stasis indices. Higher frame rates ([&ge;] 20 per cycle) offer improved quantitative accuracy, while lower temporal resolutions may remain informative for patient stratification purposes, where relative ranking is more relevant than absolute accuracy.

AO_SCPLOWBSTRACTC_SCPLOWIn this work, we publish the 3D unsteady Navier-Stokes numerical simulations and meshes of coronary arteries reconstructed from invasive X-ray coronary angiograms acquired during in the Fractional Flow Reserve versus Angiography for Multivessel Evaluation 2 (FAME 2) trial. Out of the 914 clinical images, 779 vessels are successfully reconstructed and meshed. The remaining 135 vessels have been discarded since they exhibited self-intersecting geometry during the reconstruction process. The meshes are hexahedral and all of them have the same number of vertices and identical connectivity; their high quality is demonstrated using standard mesh quality indices. The simulations are performed using the Finite Element Method (FEM) with state-of-the-art coronary boundary conditions applied at the outlet. The motivation behind this effort relies from the scarcity of publicly available numerical haemodynamics data, despite the growing interest in data-driven modeling and machine learning techniques. The database is available at the link: https://doi.org/10.7910/DVN/GPCUNS

Vector flow mapping (VFM) reconstructs left-ventricular (LV) blood velocity from color-Doppler echocardiography by combining the measured beamwise component with physical and regularizing constraints. Analysis of the discrete VFM formulation shows that the inverse problem is intrinsically ill posed: the occurrence of singular modes can be predicted from the geometry of the segmented blood-pool domain, the imposed boundary conditions, and the degree of smoothing. These modes can propagate uncertainty along entire transverse bands of the reconstructed velocity field, yet conventional VFM neither quantifies this uncertainty nor allows for correcting the blood-pool segmentation. We introduce Bayesian VFM (B--VFM), a hierarchical framework that jointly infers radial and transverse velocities, a probabilistic blood-pool mask, their spatially resolved uncertainties, and hyperparameters weighting Doppler and segmentation fidelity, mass conservation, boundary conditions, and smoothness. The discretized posterior admits a closed-form gradient and exact Hessian, enabling computationally efficient, gradient-based MAP estimation, sampling, and direct analysis of ill-posed modes. Posterior inference combines Gibbs sampling of conjugate Gamma-distributed hyperparameters with conditional maximum-a-posteriori estimation and a Laplace approximation for the high-dimensional velocity and mask fields. To accommodate systematic departures from planar mass conservation, B-VFM can learn the covariance of the planar divergence residual from an ensemble of flows and incorporate it as a structured model-discrepancy prior. Independent chains converged reproducibly, while covariance priors learned from flow ensembles illustrated how model discrepancies can be incorporated into the inference. B--VFM was evaluated using Lamb-Chaplygin dipoles under ideal conditions and with Doppler corruption, Doppler voids, and segmentation defects, and using the Hicks-Moffatt family of spherical vortices to assess violations of planar mass conservation. The method produced smooth reconstructions, localized uncertainty near unreliable measurements and regions of model inconsistency, and used flow information to correct segmentation errors. Within the tested vortex family, the data-informed planar divergence prior reduced velocity bias and mask distortion. B--VFM thus provides an uncertainty-aware reconstruction method and a flexible foundation for future VFM formulations incorporating additional priors, observations, and physical models. Future work will evaluate the method using clinical data and more complex three-dimensional benchmark flows.

Finite element (FE) head models are virtual tools to study brain biomechanics and their predictions must be numerically convergent. Previous convergence studies focused on the influence of mesh size, but the potential effect of element formulation on model convergence was often ignored. To address this, one original model with brain mesh size as 6.4 {+/-} 1.9 mm was modified to generate three derivatives with the same mesh topology but different element sizes, i.e., a coarse model (mesh size: 12.2 {+/-} 3.9 mm), a medium model (mesh size: 3.2 {+/-} 1.0 mm), and a fine model (mesh size: 1.6 {+/-} 0.5 mm). Three commonly used element formulations, i.e., reduced integration, selectively reduced (S/R) integration, and full integration, were implemented to the brain elements. These models were subjected to rotational loadings along the axial, coronal, and sagittal axes, respectively. The maximum relative displacement at representative sites and 95th percentile maximum principal strain at the whole brain level were used to evaluate mesh convergency. The results showed that the S/R integration yielded a 5% difference between the original and medium meshes, while the reduced and full integration revealed a difference over 5% even between the medium and fine meshes. This study verified that the mesh convergence of FE brain models is affected by the choice of element formulation and the S/R integration contributes to the fastest convergence behavior than the reduced and full integrations. It provided practical information on how to develop numerically convergent and computationally efficient FE brain models. HighlightsO_LIThis study verifies that the choice of element formulation affects the mesh convergence behavior of finite element brain models C_LIO_LIThis study finds the selectively reduced integration yields the fastest convergence behavior than the reduced and full integration C_LIO_LIThis study provides practical guidance on the choice of mesh density and element formulation on how to develop numerically convergent and computationally efficient finite element brain models. C_LI

Cerebrovascular pulsatility drives measurable brain tissue deformation and has been associated with ageing and a range of neurological disorders. Yet how pulsatile haemodynamic forces are transmitted through deformable cerebral arteries into the surrounding brain remains poorly understood, particularly in anatomically realistic vascular geometries. Existing computational approaches have largely treated cerebral fluid and tissue mechanics separately or relied on idealised geometries, limiting our ability to determine how vascular anatomy simultaneously governs intraluminal haemodynamics and extravascular mechanical loading. Here, we develop an image-derived three-dimensional computational framework that jointly resolves pulsatile blood flow, arterial wall deformation and surrounding brain tissue motion in representative cerebral arteries. Four arterial segments, including the middle cerebral artery, middle cerebral artery bifurcation, basilar artery and internal carotid artery, are reconstructed from high-field (5 Tesla) magnetic resonance imaging data of a healthy subject. A finite-deformation fluid-structure interaction model is established by coupling non-Newtonian blood flow, hyperelastic arterial wall and hyper-viscoelastic brain tissue. The predicted tissue response is benchmarked against in vivo magnetic resonance elastography measurements of cardiac-induced volumetric strain over a cardiac cycle. Results reveal spatially localised arterial and tissue deformation whose magnitude and distribution are strongly governed by vascular geometry and wall thickness. Among the segments examined, the internal carotid artery exhibits the largest deformation response, while reduced wall thickness increases strain transmission into the surrounding tissue. Geometrically complex regions also exhibit greater spatial heterogeneity in near-wall haemodynamic metrics. These findings demonstrate that cerebral vascular anatomy simultaneously shapes intraluminal haemodynamics and extravascular mechanical loading. By integrating image-derived vascular anatomy, coupled blood-vessel-brain mechanics and in vivo benchmarking within a unified framework, this study provides a mechanically consistent reference for healthy cerebral pulsatility and establishes a foundation for quantifying how blood-vessel-brain interactions are altered under pathological conditions.

The kidney is a critical soft-tissue organ responsible for blood filtration. Accurate constitutive models of the kidney are essential to predict tissue deformation and stress, yet existing models prescribe the strain energy function a priori, largely treat tissue variability as deterministic, and do not distinguish between cortex and medulla. Here we use Gaussian constitutive artificial neural networks to discover probabilistic strain energy functions for kidney cortex and medulla from tension, compression, and shear experiments. The medulla is approximately twice as stiff as the cortex across all three modes: The effective Young's moduli are 3.41, 5.04, and 4.48 kPa for the medulla and 1.43, 2.32, and 2.87 kPa for the cortex in tension, compression, and shear. Both regions display pronounced tension--compression asymmetry. When trained on all three modes, the network discovers strain energy functions in the second invariant alone over the tested deformation range. The cortex model selects two exponential I2 terms, while the medulla model selects an additional linear I2 term. Gaussian external weights propagate specimen-to-specimen variability into closed-form probabilistic stress predictions. Our region-specific probabilistic models capture the mechanical heterogeneity and variability of the kidney and enable more realistic simulations of kidney deformation and stress for surgical planning, needle interventions, and renal trauma.

Devireddy et al. [1] previously showed that an enthalpy-based, macroscale model simulated tissue freezing histories are in close agreement with a coupled model incorporating cellular water transport and intracellular ice formation (IIF). That comparison, however, was performed for a limited set of tissue geometries, cooling conditions, and biophysical parameters. Here, we revisit that conclusion by reimplementing the published model from its governing equations and systematically expanding the parameter space to determine when microscale processes alter the predicted macroscale freezing response. We incorporate both surface catalyzed nucleation (SCN) and volume catalyzed nucleation (VCN) of intracellular ice with cellular osmotic dehydration and examine a broad range of tissue dimensions, cooling rates, convective boundary conditions, membrane permeabilities, cell sizes, activation energies, intracellular water fractions, and nucleation rates. Linearization of the cellular water transport equation about osmotic equilibrium yields an osmotic relaxation timethat can be compared with the local residence time of tissue undergoing phase change. Their ratio defines a dimensionless group, that organizes the transition between regimes in which the uncoupled enthalpy formulation and the coupled micro and macroscale formulation produce similar or substantially different results. Across the one-dimensional parameter space examined here, the two models remain within a prescribed thermal-error bound for the dimensionless group below approximately 0.1, whereas the discrepancy increases systematically above this threshold. Model divergence also depends on intracellular ice nucleation, i.e., when nucleation is sufficiently slow, cellular water transport can become rate limiting for latent heat release. Conversely, sufficiently rapid intracellular ice formation reduces the macroscopic discrepancy by providing an additional pathway for latent-heat release. Multidimensional simulations further show that geometry alone can move a nominally conventional freezing protocol into the coupled regime. In particular, multidirectional cooling can reduce the local residence time sufficiently to produce substantial differences in predicted frozen volume even when the corresponding thermal histories appear similar. Finally, we formulate an inverse approach for estimating IIF parameters from tissue scale thermal histories. Synthetic simulations indicate that parameter recovery is feasible within a finite range of cooling conditions, providing a testable strategy for future experimental measurements. Together, these results provide a mechanistic criterion for deciding when microscale biophysical processes must be retained in macroscale tissue-freezing models.

Objective: Extracellular electric fields critically influence neural activity through both exogenous neuromodulation and endogenous ephaptic coupling. While conventional cable models efficiently simulate membrane dynamics, they fail to capture bidirectional, field-mediated interactions self-consistently, and fully coupled volumetric methods require computationally prohibitive 3D meshing. We present Cable-BEM, a hybrid wire-kernel bidomain boundary element method designed to resolve these limitations. Approach: By analytically integrating boundary integral kernels around cylindrical neuronal compartments, Cable-BEM fully couples intracellular, extracellular, and membrane dynamics while strictly retaining the highly efficient 1D degrees of freedom of traditional cable equations. The system is advanced using a semi-implicit Crank-Nicolson scheme. To overcome the dense nature of the resulting integral operators, we implement an Adaptive Cross Approximation (ACA) and Hierarchical Off-Diagonal Low-Rank (HODLR) compression scheme. Main result: The solver was rigorously validated against full-surface bidomain boundary element method (BEM) reference implementations, demonstrating tight agreement in activation thresholds (within 1.3% relative error) across diverse stimulation geometries. The ACA-HODLR compression scheme achieved substantial memory footprint reductions-by a factor of up to 4.6 for large 225-cell networks- without sacrificing numerical accuracy. Furthermore, we utilized the framework to resolve subtle, distance-dependent ephaptic interactions, successfully demonstrating the progressive phase synchronization of biophysically realistic, multi-compartment Purkinje cells. Significance: Cable-BEM provides a computationally scalable, mesh-free framework that establishes a powerful and practical foundation for investigating complex field-mediated phenomena in large-scale, multicellular neuronal networks.

Purpose: Atrial fibrosis is associated with thromboembolic risk in atrial fibrillation, but how fibrosis-related contractile impairment alters blood transport in the left atrium (LA) and left atrial appendage (LAA) remains unclear. Prior multi-physics simulations coupled electromechanical (EM) and computational fluid dynamics (CFD) models to show that fibrosis impairs LA function and flow, but were limited to booster-only contraction and could not separate the effects of fibrosis burden from its spatial distribution. Here we examined how fibrosis burden and spatial pattern independently relate to long-term LA and LAA hemodynamics in full-cycle, multi-beat simulations on a common LA anatomy. Methods: Three fibrosis patterns derived from late gadolinium enhancement MRI were considered: one from the patient providing the common anatomy and two from additional patients, with the latter mapped onto the common anatomy via universal atrial coordinates. All three patterns were applied globally at matched burdens of 15.6% and 31.2%, yielding six global-fibrosis cases and a non-fibrotic reference. Two patterns were also confined to the LAA at three local burdens (10%, 25%, and 50%), yielding six LAA-fibrosis cases. A 3D EM model for the LA anatomy was coupled to a 0D circulatory model, so that the simulations resolved the reservoir, conduit, and booster phases. Full-cycle EM-derived wall motion drove subsequent multi-beat CFD simulations. Region-specific emptying fraction (EF), non-dimensional kinetic energy (KE), and the 90th percentile of blood residence time (TR) were evaluated for the LA and LAA. Results: Left ventricular stroke volume remained nearly constant across cases, while increasing global fibrosis burden progressively reduced LA and LAA EF (LA: 36% to 30-31%; LAA: 44% to 28-38%) and monotonically increased LA TR across all patterns. Non-dimensional LA KE showed an approximately quadratic association with LA EF across patterns, whereas LAA KE was more scattered but still increased with LAA EF. In contrast, the relationship between the 90th percentile of LAA TR and LAA EF showed pattern-specific, non-monotonic trends. Conclusions: Fibrosis burden and spatial distribution jointly shape LA and LAA hemodynamics, but through different mechanisms. LAA EF is a robust functional descriptor of chamber-averaged kinetic energy regardless of fibrosis pattern, whereas LAA blood stasis depends additionally on the spatial distribution of fibrosis beyond what bulk emptying fraction predicts. These findings support incorporating imaging-derived fibrosis pattern, and not only burden, alongside functional metrics when assessing thrombogenic risk in the LAA.

Healing of the gastrointestinal (GI) tract following surgical interventions is governed by coupled mechanobiological processes and may lead to severe complications such as perforation or stricture formation. A first comprehensive mechanobiological modeling framework of GI healing is herein proposed. The framework addresses two clinically relevant applications of bioprinting: endoscopic treatment of GI tissue defects and bioprinting-assisted surgical anastomosis. Tissue remodeling is based upon a hyperelastic material model with fiber dispersion coupling the spatiotemporal reaction diffusion dynamics of fibroblasts and cytokines with collagen deposition, reorientation, and plastic deformation. An active stress formulation drives fibroblasts and myofibroblast contraction while a multiplicative decomposition approach is adopted for collagen remodeling as a permanent change. The governing equations are discretised in time with an implicit Euler scheme, finite element in space and solved by a staggered algorithm implemented in the open-source platform FEniCS. Numerical simulations show that the initial collagen fibers content and orientation strongly influence fibroblast infiltration and scar contracture formation. The study quantifies the advantage of having a collagen bioprinted architecture aligned with tissue fibers, to reduce structural remodeling distortion. The model also predicts a strong correlation between collagen deposition and luminal burst pressure during the recovery of mechanical integrity in GI anastomosis.

Reduced-order models are increasingly used to study cardiac physiology and inform patient-specific therapies. However, a model's prediction is only as reliable as its underlying parameters: representative model parameterization is essential to reflect the physiology of the populations these models are meant to represent, including biological sex. Most current models are parameterized from male or sex-agnostic data and/or focus on specific pathologies. Therefore, the goal of this work is to establish a formal parameter estimation pipeline for deriving reduced-order cardiovascular model parameter ranges that are physiologically representative of healthy women and men. We calibrated a closed-loop reduced order model of the heart and circulation separately for healthy female and male populations, using data pooled from eleven healthy cohorts. To account for parameter sensitivity and identifiability, we employed a three-stage parameter subset reduction pipeline: global sensitivity analysis (Sobol's method), collinearity screening (Fisher information matrix), and profile-likelihood identifiability analysis. Sex-specific distributions of parameters that were deemed sensitive and identifiable for each sex, ten for women and 9 for men, were obtained by Hamiltonian Monte Carlo. All the calibrated parameters showed less than 80% overlap between sexes, with the smallest overlap observed in some of the most influential parameters, such as stressed blood volume. Comparing simulations of the calibrated models against allometrically size-matched simulations showed that body size explained some, but not all, of the sex differences. The resulting parameter distributions provide reference ranges usable in future mechanistic and patient-specific simulations to contribute to more inclusive cardiovascular modeling.

Pulmonary hypertension (PH) is characterized by progressive structural and mechanical remodeling of the pulmonary vasculature, yet few computational frameworks directly link disease mechanisms to longitudinal progression and therapeutic response. In this study, we utilized a multiscale pulmonary arterial growth and remodeling (G&R) framework to capture evolving functional metrics from rat models of PH. This framework couples morphometric tree hemodynamics, constrained mixture theory-based wall mechanics, and maladaptive cellular remodeling. Disease progression was driven by three mechanistically interpretable parameters governing excess smooth muscle production, remodeling activation, and passive stiffening. These parameters were calibrated to longitudinal monocrotaline (MCT) measurements of pressure, wall thickness, and stiffness from prior work using a multiobjective optimization. To show the predictive value of this model, we simulated therapeutic intervention within the same disease-specific framework by using functional cell-level responses to therapy to inform changes in parameter values. Calibration to the study-specific MCT dataset reproduced the temporal increases in pressure, wall thickness, and stiffness, demonstrating that the model could capture multiple features of vascular remodeling simultaneously, with R2 values of 0.81, 0.83, and 0.95, respectively. Simulated treatment reduced pressure, wall thickness, and stiffness. Predicted pressure and wall-thickness responses agreed closely with the corresponding experimental treatment effects, whereas stiffness recovery was overpredicted, suggesting that additional mechanisms may contribute to persistent vascular stiffening after intervention. The framework also captured the overall progression of pulmonary pressure increases across both aggregated MCT and Sugen-hypoxia datasets, suggesting utility across studies and animal models. This work outlines a physics-based, multiscale framework that simulated quantities of direct clinical interest in a mechanistically interpretable platform for linking pulmonary vascular remodeling and treatment response. It supports comparisons across experimental phenotypes and interventions while identifying where constitutive refinements are needed to improve predictive capability across phenotypes.

This study presents a physiologically detailed biomechanical model of the mouse distal forelimb that incorporates intrinsic musculature, tendon routing, and digit-level skeletal anatomy, features simplified or omitted in existing musculoskeletal models. Using high-resolution anatomical reconstruction and computational modeling, we created a physiological representation of the wrist and digits capable of simulating complex forelimb movements. The model enables simulation of coordinated distal forelimb movement and digit-level muscle behavior during grasping-related tasks. Simulations were performed for multiple tasks, including grasping, grasping with supination, wrist flexion, and digit I flexion, with analysis focused on the grasping task due to its integration of both intrinsic and extrinsic musculature. Model performance was evaluated through comparisons of marker trajectories between torque-driven reference motion and muscle-driven simulations, temporal shuffle control, and comparisons between experimentally recorded electromyography (EMG) activity and model-predicted muscle excitation profiles. The model successfully reproduced coordinated distal forelimb kinematics, demonstrated strong agreement between torque-driven and muscle-driven simulation approaches, and generated physiologically plausible muscle excitation patterns consistent with experimentally observed EMG activity during grasping-related movement. These findings establish the model as a framework for studying fine motor control, neuromuscular coordination, and movement-related impairments in mice while providing a foundation for future investigation of neurological disorders and their underlying biomechanical mechanisms.

A densely vascularized tissue can still be starved of oxygen. Tissue oxygenation depends heavily on vessel arrangement, blood flow distributions, and cellular metabolism. Resolving these effects from individual constituent vessels to whole organs has historically been computationally prohibitive. We introduce the Cardiovascular Analysis and Simulation of Circulatory Advection, Diffusion, and Exchange (CASCADE) toolkit, which makes vessel-resolved oxygen transport tractable at the organ scale through a novel screened Green's-function formulation. On laptop-class hardware, CASCADE simulates networks containing thousands of vessels in under one second and a full biventricular human heart model in about two minutes. We validate spatial oxygen and cell-survival predictions against experiments in perfused engineered tissues and matched finite-element simulations. In engineered tissues, sweeps across thousands of network designs and physiological conditions reveal when vessel density, oxygen supply, or network architecture limits viability, establishing design rules for sustaining thick constructs. In a human heart model, CASCADE translates coronary occlusion, anemia, polycythemia, and hypoxemia into distinct myocardial oxygenation patterns, revealing where oxygen delivery fails and how much tissue is at risk.

Cerebral aneurysms are abnormal outpouchings of arteries within the brain and occur in [~]1 in 30 adults. Their initiation, growth, and rupture have been linked to focal blood flow abnormalities--often termed "disturbed" or "hostile" hemodynamics--but commonly-used hemodynamic metrics yield conflicting associations with pathology and lack a unifying mechanistic interpretation. Building on a theoretically-grounded link between wall shear stress and near-wall vorticity, we hypothesized that a topology-based description of near-wall flow can operationalize the concept of hostile hemodynamics in a reproducible way. Inspired by atmospheric tornadic phenomena, we sought a principled taxonomy of coherent near-wall fluid structures with potential mechanobiological and clinical implications. Using high-fidelity computational fluid dynamics simulations in anatomically realistic geometries, we identified coherent near-wall fluid structures whose organization mirrors well-studied atmospheric phenomena: tornado-like columnar rotating cores; downburst-like nonrotating wall-impinging jets with tangential outflow, roll-cloud-like tangential vortices; and mixed configurations. These tornadic events on the aneurysm luminal surface were identified from wall shear stress topology, consistent with its theoretical connection to near-wall vorticity kinematics. The presence of tornadic phenomena--and their imprints on the aneurysm wall--was independently observed in vivo using 4D flow magnetic resonance imaging. By translating concepts from atmospheric physics into vascular biomechanics, this topology-based framework yields a unified mechanistic language for describing near-wall hemodynamics, resolving blood flow complexity into interpretable and reproducible coherent fluid structures, enabling standardized hemodynamic phenotyping, and supporting hypothesis-driven studies of aneurysms and other cardiovascular diseases where greater fluid-mechanical specificity and interpretability may strengthen links between mechanobiology and clinical risk.

Cancer survival rates increase with earlier diagnosis. However, many cancers are diagnosed only after symptoms develop, in the later stages of cancer progression. We need a diagnostic tool that can detect cancer earlier. Circulating tumour cells (CTCs) are cells that detach from the main tumour and can enter the bloodstream. The capture and analysis of CTCs can provide an early indicator of cancer. However, rarity of CTCs in blood makes their efficient capture challenging. Inertial microfluidics can be utilised for a label-free separation of CTCs from white blood cells (WBCs) by manipulating cells in the microchannel based on cell properties, achieving good separation performance. Inertial microfluidic devices mainly rely on size-based separation. However, residual WBC carryover limits complete separation. Here, we show that cell mechanical heterogeneity plays a large role in the separation process and should be considered in the design of these devices. We simulate the migration dynamics of CTCs and WBCs in a straight microchannel using a 3D lattice-Boltzmann-immersed-boundary-finite-element solver. Our results demonstrate that WBC migration behaviour changes depending on the deformability of a CTC. The size and deformability of the cells have been shown to determine the cell-specific heteroclinic orbits leading to different interaction types, these interactions alter WBC migration resulting in more/less WBC carryover. This work highlights that both cell-specific heteroclinic orbits and single cell migration rates should be considered in the design of inertial microfluidics for separation.

Muscle architecture shapes muscle function and changes with age, growth, training and disease, yet quantifying three-dimensional (3D) muscle architecture in vivo remains challenging. We introduce a hybrid fascicle tractography approach for freehand 3D ultrasound data that accurately reconstructs 3D muscle fascicles with respect to an objective, anatomically relevant coordinate system defined by the muscle's central aponeurosis. The hybrid approach combines Hessian-based fascicle detection with wavelet-based refinement to generate volumetric fascicle orientations. In a synthetic dataset with known ground truth, fascicle orientations and lengths were estimated with errors of [&le;]2{degrees} and ~1.5%, respectively. In vivo, the approach detected physiologically plausible fascicle lengthening in the human tibialis anterior following a passive plantar flexion rotation, whereas diffusion tensor imaging of the same muscle did not. The proposed method enables anatomically relevant, objective and non-invasive quantification of 3D muscle architecture in vivo, providing a practical framework for applications in clinical and applied muscle physiology.